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£SJ ■ Abstract 

Using coherent-state techniques, we prove a sampling theorem for Majorana's (holomor- 



phic) functions on the Riemann sphere and we provide an exact reconstruction formula 
as a convolution product of N samples and a given reconstruction kernel (a sinc-type 
function). We also discuss the effect of over- and under-sampling. Sample points are 
roots of unity, a fact which allows explicit inversion formulas for resolution and overlap- 
ping kernel operators through the theory of Circulant Matrices and Rectangular Fourier 
Matrices. The case of band-limited functions on the Riemann sphere, with spins up to 
J, is also considered. The connection with the standard Euler angle picture, in terms of 
spherical harmonics, is established through a discrete Bargmann transform. 

X. 

c3 ■ 1 Introduction 

Despite the fact that the Fast Fourier Transform (FFT), in the setting of Abelian harmonic 
analysis (i.e. the well known Cooley-Tukey algorithm pQ for time series analysis), has been 
extensively studied in both the theoretical and applied literature, the number of algorithms for 
the efficient computation of Fourier transforms associated with non-Abelian groups and their 
homogeneous spaces is still relatively scarce. 

For finite non-Abelian groups, like the symmetric group S n , the reference [2] provides effi- 
cient algorithms to compute Fourier transforms. For the two-dimensional sphere § 2 , the efficient 
computation of Fourier transforms of band-limited functions (those functions in L 2 (S 2 ) which 
expansion requires only spherical harmonics of angular momentum at most J) has been achieved 
in, see for instance, Refs. [3j IU [U EJ . In reference [3], the authors develop a sampling theo- 
rem on the sphere, which reduces the computation of Fourier transforms and convolutions of 
band-limited functions to discrete (finite) calculations. Here, band-limited functions on S 2 , of 
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bandwidth J, are expanded in terms of spherical harmonics and sampled at an equiangular grid 
of 4J 2 points. 

The point of view followed in these references is a group theoretic one. In this setting, 
the FFT on S 2 is an algorithm for the efficient expansion of a function defined on the sphere 
§ 2 = SO{3) / SO{2) in terms of a set of irreducible matrix coefficients for the special ortogonal 
group in three dimensions, G = 50(3), which, in this case, are the standard family of spherical 
harmonics. 

In this article we consider the group G = SU{2) (double cover of SO (3)), which allows for 
(extra) half-integer angular momenta (spin). Moreover, we shall work in a different (holomor- 
phic) picture and use, instead of spherical harmonics, other system of (less standard) orthogonal 
polynomials: "Majorana's (holomorphic) functions" 0[8] on the Riemann sphere C = CU {00} 
(one-point compactification of the complex plane). We shall give a generalized Bargmann Trans- 
form [9] relating both pictures (representations): the "holomorphic" one and the "standard" one, 
which is a particular case of coherent-state transform |10[ [TTJ . Then we shall choose in C the 
roots of unity as sampling points, so that the sampling of the coherent-state overlap (or Repro- 
ducing Kernel) has a "circulant" structure [12] . Using the properties of the Rectangular Fourier 
Matrices (RFM) and the theory of Circulant Matrices we will be able to invert the (sampled) re- 
producing kernel B and provide a reconstruction formula for Majorana's (holomorphic) functions 
on the Riemann sphere. The inversion formula is accomplished through an eigen-decomposition 
B = J-DJ 7 ^ 1 of £>, where T turns out to be the standard discrete fourier transform matrix. This 
fact allows for an straightforward fast extension of the reconstruction algorithm. 

In order to keep the article as self-contained as possible, we shall introduce in the next two 
sections general definitions and results about coherent states and frames based on a group G, and 
the standard construction of spin coherent states for the case G = SU (2). We address the reader 
to Refs. [101 [TTj [TJJ [15] for more information. In Section 0] we provide sampling theorems and 
reconstruction formulas for Majorana's functions on the Riemann sphere, discussing the effect of 
over- and under-sampling and the analogies with the so called "covariant interpolation" . We also 
discuss the case of band-limited functions. In Section [5] we provide explicit expressions (discrete 
Bargmann transforms) which connect our "complex holomorphic picture" and the standard 
"Euler angle picture", and we discuss some obstructions that arise. Appendixes A and B are 
devoted to a brief review on rectangular Fourier matrices and circulant matrices, respectively. 

2 A brief on Coherent States and Frames 

Let us consider an unitary representation U of a Lie group G on a Hilbert space (T~t, (■{•})■ 
Consider also the space L 2 (G,dg) of square-integrable complex functions \P on G, where dg = 
d(g'g), V</ G G, stands for the (let us say) left-invariant Haar measure, which defines the scalar 
product: 



Jg 

A non-zero function 7 G 7i is called admissible (or fiducial vector) if T(g) = (U(g)j\^) G 
L 2 (G, dg), that is, if 




(1) 
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Let us assume that the representation U is irreducible, and that there exists a function 7 
admissible, then a system of coherent states (CS) of H associated to (or indexed by) G are 
defined as the functions in the orbit of 7 under G: 

l 9 = U(g) 1 , g£G. (3) 

We can also restrict ourselves to a suitable homogeneous space Q = G/H, for some closed 
subgroup H. Then, the non-zero function 7 is said to be admissible mod(-ff, a) (with a : Q — > G 
a Borel section), and the representation U square integrable mod(-ff, a), if the following condition 
holds: 

r \(U(a(q)) 1 \4,)\ 2 dq<cx>, V^«, (4) 



where dq is a measure on Q "projected" from the left-invariant measure dg on the whole G. 
The coherent states indexed by Q are defined as J a (q) = U{a{q))^, q 6 Q, and they form an 
overcomplete set in TL. 

The condition dH) could also be written as an "expectation value" 

0< / \{U{a{q)) 1 \ij)\ 2 dq=^\A u \^) < 00, V^W, (5) 
JQ 

where A a = Jq \j a (q)) (la(q)\dq is a positive, bounded, invertible operator[j] If the operator A^ 1 
is also bounded, then the set S a = {|7o-( g ))) <? £ Q} is called a frame, and a izg/zi frame if ^ is 
a positive multiple of the identity, A a = XI, A > 0. 

To avoid domain problems in the following, let us assume that 7 generates a frame (i.e. that 
A" 1 is bounded). The CS map is defined as the linear map 



T 7 : H — * L\Q,dq) 



(6) 

Its range L^(Q, dq) = T 7 (7i) is complete with respect to the scalar product ($|^) 7 = (<£|T 7J 4 a : 1 T~ ll I') 

and T 7 is unitary from H onto L 2 (Q,dq). Thus, the inverse map T" 1 yields the reconstruction 
formula: 

i> = T- 1 %= f y^A-^^dq, ^> y GL 2 (Q,dq), (7) 
J Q 

which expands the signal ip in terms of CS A~ 1 ^ a ^ q \ ) with wavelet coefficients ^^{q) = [T^tp](q). 
These formulas acquire a simpler form when A a is a multiple of the identity, as is for the case 
considered in this article. 

When it comes to numerical calculations, the integral A a = jn\la(q)){lcj(q)\dq has to be 
discretized, which means to restrict oneself to a discrete subset Q C Q. The question is whether 
this restriction will imply a loss of information, that is, whether the set S = {\qk) = \lcr(q k ))-, Qk £ 
Q} constitutes a discrete frame itself, with resolution operator: 

A=J2\lk)(lk\- (8) 



^In this paper we shall extensively use the Dirac notation in terms of "bra" and "kets" (see e.g. [161 115] ). 
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The operator A need not coincide with the original A a . In fact, a continuous tight frame might 
contain discrete non-tight frames, as happens in our case (see later on Sec. 0|). 

Let us assume that S generates a discrete frame, that is, there are two positive constants 
< b < B < oo (frame bounds) such that the admissibility condition: 

&||VH 2 <I (qk\1>)\ 2 < BM\\ 2 (9) 

9fcGQ 

holds Vip £ TC. To discuss the properties of a frame, it is convenient to define the frame (or 
sampling) operator T : TC — > I 2 given by T(ip) = {{q k \^), q k £ Q}- Then we can write A = T^T, 
and the admissibility condition (|9|) now adopts the form: 

bl < T ] T < BI, (10) 

where I denotes the identity operator in TC. This implies that A is invertible. If we define the 
dual frame {\q) = A~ 1 \q)}, one can easily prove that the expansion (reconstruction formula): 

m = E famm (n) 

gk&Q 

converges strongly in TC, that is, the expression 

T+T = E \Qk)(q k \ = T\T^ = ^ \lk}(q k \ = I (12) 

provides a resolution of the identity, where T ( + = (T^T)~ l T^ is the (left) pseudo-inverse (see, for 
instance, [13]) of T (see e.g. [HI [15] for a proof, where they introduce the dual frame operator 
f = (7J+) 1 " instead). 

It is interesting to note that the operator P = TT t + acting on £ 2 is an orthogonal projector, 
as can be readily verified, onto the range of T. 

We shall also be interested in cases where there are not enough points to completely recon- 
struct the signal, i.e. under sampling, but a partial reconstruction is still possible. In these cases 
S does not generate a discrete frame, and the resolution operator A would not be invertible. 
But we can construct another operator from T, B = TT^ , acting on £ 2 . 

The matrix elements of B are Bki = (<7fc|<ft), therefore B is the discrete reproducing kernel 
operator, see eq. (|24p . If the set S is linearly independent, the operator B will be invertible 
and a (right) pseudo- inverse can be constructed for T, = (TT^)^ 1 , in such a way that 
= Ip.. As in the previous case there is another operator, P$ = T+T acting on TC which is 
an orthogonal projector onto the subspace spanned by S. A pseudo-dual frame can be defined 
as 

\q k )= Y,iB~ l )lk\m) (13) 

qieQ 

providing a resolution of the projector P$ : 

T+T = £ \q k )(q k \ = Tt(T+)t = £ \q k )(q k \ = P s (14) 
g^eS q k eQ 

Using this, an "alias" of the signal is obtained, 

gfceQ 
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which is the orthogonal projection of \ip) onto the subspace spanned by S, \ip) = Ps\ift}- An 
example of this can be found in Sec. 14.1.21 

The two operators A and B are intertwined by the frame operator T, TA = BT. If T is 
invertible, then both A and B are invertible and = = T~ 1 . This case corresponds to 
critical sampling, where both operators A and B can be used to fully reconstruct the signal. 

3 Representations of SU(2): Spin Coherent States 

Spin coherent states are based on the familiar Lie algebra of self-adjoint spin operators su(2) = 
{Ji, J2, J3} satisfying the standard angular momentum commutation relations [J±, J%] = 1J3 
(plus cyclic permutations). The spin operators Jj~ are the infinitesimal generators of the group 
SU(2) and any element U £ SU (2) can be represented in any of several parameterizations: 

Euler angles : U(9, </>, ip) = e^^e^^e'^ 3 , 
Canonical coordinates : U{a u a 2 ,a 3 ) = e i{ai Jl+a2 J ^+ a ^\ (16) 
Complex Coordinates : U(z, z, ip) = J\fe zJ ~ e~ zJ+ e~ lipj3 , 

where J± = J\ ±i J 2 are the raising and lowering ladder operators and M a suitable normalization 
factor that makes U unitary (its explicit expression will be given later). The complex ("Hopf 
fibration") coordinates z,z are related to the stereographic projection of the 2-sphere on the 
complex plane, z = tan(|)e i<?i , and this is usually referred to as the Riemann sphere. 

Unitary irreducible representations of the Lie algebra su{2) are (2s + l)-dimensional, where 
s = 0,l/2, 1,3/2,... is a half- integer parameter (spin or angular momentum) that labels each 
representation. Each carrier space Tt s ~ C 2s+1 is spanned by the common angular momentum 
orthonormal basis B(7i s ) = {\s, m), m = —s, . . . , s} (in bra-ket notation) of eigenvectors of J3 
and the Casimir operator J 2 = j\ + jf + Jf. The standard action of the spin operators Jk on 
this basis vectors is: 

Js\s,m) = m\s,m) , J±\s,m) = \/ (s =F m)(s ± m + l)|s, m ± 1). (17) 

Note that the subgroup U(l) C SU(2), generated by J3, stabilizes any basis vector up to 
an overall multiplicative phase factor (a character of U(l)), i.e. e- ivJ3 |s,m) = e~ imip \s,m). 
Thus, letting Q = SU(2)/U(l) = S 2 and taking the Borel section a : Q -> G with a(<j>,9) = 
(9, (f>,cp = 0),OT a(z, z) = (z, z, 0), we shall define, from now on, families of covariant (see pj)]) CS 
mod(t r (1), a). In simple words, we shall set if = and drop it from the vectors: U (9, (ft, <p)\s, m) 
and U (z, z, <p)\s, m). 

Therefore, we have different characterizations of spin coherent states according to distinct 
choices of parameterizations. We shall concentrate in the complex and Euler angle parameteri- 
zations. 

3.1 Euler Angle Characterization: Spherical Harmonics 

For any choice of fiducial vector \y) = \s,m) the set of coherent states \9,(ft;m) = U(9,(j))\^) is 
overcomplete (for any m) in TC S . They can be easily computed by exponentiating the relations 
(|17p . This set of coherent states is also a tight frame with 

J= 2s + i r i0 ;m \/0 ^i^Q n 8 ) 

4vr 7§2 
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a resolution of unity and dfl = sin 9d0d<fi the standard invariant measure on the 2-sphere. 

For the particular case of integer spin s = j and fiducial vector m = 0, the standard spherical 
harmonics YJ™(#,(/j) arise as the irreducible matrix coefficients: 



(9,<f>;0\j,m) = (jAU{0,ct>)- l \j,m) = ^1^(0,0), (19) 

or, in other words, the components of spin coherent states \0, (f>;0) over the orthonormalized 
basis {\j, m)}. Thus, for a general angular momentum j state we have the standard spherical 
harmonic decomposition [the wavelet coefficients ©]: 



*(e, </>) = (9, </>; 0|V) = J £ ^ m Y™{6, <j>), (20) 

m=-jf 

with Fourier coefficients Vv« = {j,m\ip). 

3.2 Complex Holomorphic Characterization: Majorana functions 

In this case we shall use I7) = \s,m = s) as fiducial vector (i.e. the highest weight vector), so 
that J+I7) = and the coherent states \z) = U(z,z)\ / y) = Ne zJ ~\^) are holomorphic (only a 
function of z), apart from the normalization factor M, which can be easily computed by using 
the relations (fT7|) giving J\f = (1 + zz)~ s . The frame {\z), z £ C} is also tight in Tt s , with 

2s + 1 /"..., d 2 z . , 

kX^l Ti 1 -N2 » ( 21 ) 



vr i c IMI (l + ^)2' 

where we denote d 2 z = dRe(z)dlm(z). 

For this case, the decomposition of de CS \z) over the orthonormalized basis {|s,m)} gives 
the irreducible matrix coefficients: 

(z\s,m) = ( S , S \U(z,z)- 1 \s,m) = ( s l s m ) 1/2 (l + zzy s z s+m 

= (l + z-zr s r™(z), (22) 

where now T™(z) is just a monomial in z times a numeric (binomial) factor. A general spin 
s state is represented in the present complex characterization by the so called Majorana 
function [3 [8] : 

s 

9(Z) = (Zty) = (1 + ZZ)- S J2 ^rn??(z) = Aff(z), (23) 

m=—s 

which is an anti-holomorphic function of z (in this polynomial) H 

Note that the set of CS {\z)} is not orthogonal. The CS overlap (or Reproducing Kernel) 
turns out to be: 

Cl + z 'z) 2s 

C(z, z>) = (z\z>) = (24) 
(1 + zz) s {l + z'z'Y 



'''Here we abuse of notation when representing the non-analytic function ^(z, z) simply as 9(z), which is indeed 
anti-holomorphic up to the normalizing, non-analytic (real), pre-factor TV = (1 + zz) ~ s . Usually, this pre-factor 



is absorbed in the integration measure in (|21|l . If we choose the lowest weight fiducial vector I7) = \s,m 
we would obtain proper holomorphic functions f(z). 
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This quantity will be essential in our sampling procedure on the Riemann sphere. 

For completeness, let us provide an expression which interplays between both characteriza- 
tions of CS for integer spin s = j. It is given by the CS (or Bargmann-like) transform (see e.g. 

mm- 

K(B,hz) = (6A\z)= E (0A\j,m)(j,m\z) = (l + zz)-^J^^ £ YJ" (0 , <j>)T? (z) 

m=—j m=—j 

= (1 + zz)- j ^^A^i n B e -^ + 2zcos6 - z 2 sm6e^) j , (25) 

which can be seen as a generating function for the spherical functions YJ a (6,(p) when we drop 
the normalization factor N from the last expression. 

4 Sampling Theorem and DFT on § 2 

Sampling techniques consist in the evaluation of a continuous function ( "signal" ) on a discrete set 
of points and later recovering the original signal without loosing information in the process, and 
the criteria to that effect are given by various forms of Shannon's Sampling Theorem. Basically, 
the density of sampling points must be high enough to ensure the reconstruction of the function 
in arbitrary points. We shall concentrate on fixed spin holomorphic (Majorana's) functions and 
sample them at the roots of unity. 

4.1 Single spin case 

Let us first restrict ourselves to functions in 7i s , i.e. with well-defined spin or angular momentum 
s. In this case there is a convenient way to select the sampling points in such a way that the 
resolution operator A and/or the reproducing kernel operator B are invertible and explicit 
formulas for their inverses are available. These are given by the N roots of unity in the complex 
plane, N £ N, which would be associated, by inverse stereographic projection, to a uniformly 
distributed set of points in the equator of the Riemann sphere. The choice of roots of unity 
is made for convenience, since the ./V roots of any non-zero complex number would also be 
valid, and would correspond to different parallels in the Riemann sphere, but then the formulas 
obtained are less symmetrical than the ones corresponding to roots of unity. The most important 
reason to select roots of unity is that they are associated with the discrete cyclic subgroup 
C U(l) C SU (2). The choice of N = 2s + 1 roots of unity corresponds to a critical sampling. 
We shall also discuss the consequences of over-sampling, with N > 2s + 1, and under- sampling, 
with N < 2s + 1, in the following subsections. 

4.1.1 Over-sampling and critical sampling 

In the case of over-sampling the set S generates 7i s , and the resolution operator A = T^T is 
invertible. The case of critical sampling is a particular case of this and therefore the following 
discussion also applies to it. 

The previous statements are formalized by the following: 
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Lemma 4.1. Let Q = {z k = e 2nik / N , N > 2s + 1, k = 0, . . . , N - 1} be the discrete subset of 
the homogeneous space Q = SU(2)/U{1) = S 2 = C made of the N roots of unity. The discrete 
set of CS S = £ Q} constitutes a discrete frame in 7i s and the expression 



N-l N-l 



hs+i = \zk){zk\ = \zk){zk\ (26) 



fc=o k=0 



provides a resolution of the identity in H. s . Here \zk) = A 1 \zi c ) , k = 0, . . . , N — 1 , denotes the 
dual frame, and the resolution operator, A, is diagonal in the angular momentum orthornormal 
basis B{H S ), A = diag(A , . . . , A 2s ), with X n = ( 2 *) , n = 0, . . . ,2s. 

Proof. First, from eq. (|22p the expression for the matrix elements of T can be obtained, 

%n = (z k \s,n-s) = 2- s ^J]^)e~ i ^ , k = 0, 1,. . . ,N-1, n = 0, 1, . . . , 2s. Then, the resolution 
operator turns out to be: 

N-l , N-l 



A nm = ^C4 n )*T fem = 2- 2 y e 2 ^"^ N = N2- 2s ( 2 :)8 nm , (27) 

fc=0 fc=0 
where we have used the well known orthogonality relation: 

e (e 2 ~r - { S;*:;:z d N } - »*- ■ ^ 

k=0 k J 

since N > 2s + 1. Therefore ^4 is diagonal with non-zero diagonal elements, thus it is invertible 
and a dual frame and a (left) pseudo-inverse for T can be constructed, = A~ l T\ providing, 
according to eq. (fl2j) . a resolution of the identity. ■ 

Remark 4.2. It is interesting to rewrite this proof in terms of Rectangular Fourier Matrices 
(see Appendix |A|) . Let D = diag(Ao, • • • , X 2s ) be a diagonal (2s + 1) x (2s + 1) matrix, then 
T = Fn2s+1 D l l 2 = o ln,2s+i ° D 1 / 2 . From this the expression of A = T^T = D is readily 



D 












(see 

NxN 



recovered, and also B = TT^ is seen to be B = F N D^ where £> T 

Appendix [A]). Note that B is a singular NxN matrix with only 2s + 1 non-zero eigenvalues 

rt 

N 



A n , n = 0, 1, . . . , 2s, and that they coincide with those of A. In fact, B = A^ J^tr. 



Lemma 4.3. Under the conditions of the previous lemma, the operator P = TT, + = Fn Pis+i 
is an orthogonal projector onto a (2s + 1)- dimensional subspace of C N , the range of T . 
Proof. By direct computation (and using Appendix A), 

P = TT+ = TA- l T^ = r N i N}2s +i D 1 ' 2 D- 1 D 1 ' 2 p 2s+lyN F ] N = T N l n , 2s+ i P2s+i,n f N 
= FnP2 S +i^ n , (29) 



where P 2s +i = {hs+i) ] 



hs+i 












projector, unitarily equivalent to P 2s +\ and that PT = T. 



This clearly shows that P is an orthogonal 

NxN 
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Theorem 4.4. (Reconstruction formula) Any function if> S TC S can be reconstructed from 
N > 2s + 1 of its samples (the data) = (zj-\if>), at the sampling points Zk = e 2mk / N , k = 

0, . . . , N — 1, by means of: 

N-l 

*(z) = <*#> = Y, ^(^S^fc 1 ), (30) 

k=0 

where 

2 s 1 - z 2s+l 

H( a r) = -(l + ^)- n - r - (31) 

p/ays i/ie ro/e of a "sinc-type function". 

Proof. From the resolution of the identity (|26p . any ?/> G TC S can be written as \if)) = 
^2k=o ^(zk)\zk), and therefore W(z) = (z\ip) = ^2^=0 ^{ z k){z\h)- Using that \z k ) = A'^Zk), 
we derive that 

w**) = ~m £ A « 1/2e2mfcn/A ^i s < n ~ s ) = £ w 1 )" = 3 w 1 ) > * = o, i, . . . , jv-i , 

V iV n=0 n=0 

(32) 

where eq. (l22j) has been used. ■ 

Remark 4.5. It is interesting to note that eq. (|30p can be interpreted as a Lagrange-type 
interpolation formula, where the role of Lagrange polynomials are played by the functions 
Lk(z) = ^(zz^ ), satisfying the "orthogonality relations" Lf.(zi) = E(ziz^ ) = Pik, where P 
is the projector of Lemma 14.31 In the case of critical sampling, N = 2s + 1, the usual result 
L>k{zi) = dik is recovered, but for the strict oversampling case, ./V > 2s+l, a projector is obtained 
to account for the fact that an arbitrary set of overcomplete data ^(z^), k = 0, . . . , N — 1, can 
be incompatible with \if>) &H S . 

A reconstruction in terms of the Fourier coefficients can be directly obtained by means of 
the (left) pseudo-inverse of the frame operator T: 

Corolary 4.6. The Fourier coefficients a m of the expansion = Ylm=-s am \ s i m ) °f anv 
if) € TC S in the angular momentum orthonormal basis B(7i s ) can be determined in terms of the 
the data ^(z^) = {zkW) o,s: 

ns 

a n - s = ^( 2 n r 1/2 £ nzk)e 2 * ikn/N ,n = 0,...,2s. (33) 

fc=0 

Proof. Taking the scalar product with (zj.\ in the expression of \ip), we arrive at the over- 
determined system of equations 

2s 

E^fcn a n-s = *(-Zfc), %n = (Zk\s,n- s), (34) 
n=0 

which can be solved by left multiplying it by the (left) pseudo-inverse of T, = (T^T) _1 T^ = 
A~ X T^ . Using the expressions of A -1 = diag(Ap , A^j -1 , . . . , A^ 1 ), given in Lemma I3~TI and the 
matrix elements 7fc n , given by the formula (|22p . we arrive at the desired result. ■ 
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Remark 4.7. Actually, using vector notation, we have Ta = 9, where a = (a_ s , . . . ,a s ), and 
\& denotes the vector of samples ^(zk) , k = 0, . . . , N — 1. Using the (left) pseudo-inverse of 
T we can solve it obtaining a = D~ 1 / 2 p2 S +i,N Fn'® , which coincides with eq. ([33]) . Note also 
that the last expression is a map from to C 2s+1 ~ TL S due to the presence of the projector 
P2s+i,n (see Appendix [A]), and this prevents the appearance of infinities in the reciprocal of the 

binomial coefficient f^f) with n > 2s. This is clearer if we apply T to the expression of a 

to obtain Ta = Tn ln,2s+i D 1 / 2 D^ 1 / 2 P2 S +i,n = P^, that is, the data ^ should be first 
projected in order to obtain a compatible set of data. 
Next we provide an interesting expression. 

Proposition 4.8. If we define the "dual data" as T(k) = (zk\ip), then they are related to the 
data fy(k) = ^>(zk) = {zk\ip) through the convolution product: 

N-l 



T(k) = [A **](£) = A(fc-/)*(0, (35) 

1=0 

where A(k) (the filter) turns out to be the Rectangular Fourier Transform of 5 = (An , . . . , A^ 1 ), 
i.e.: 

A(k) = [T N , 2s+1 5](k) = J= £ = i E ffl-'e-*™*^. (36) 

V-'" „_n n 



n=0 n=0 

Proof. Applying ([26j) to tp we obtain: 

JV-l 



fc=0 



Taking the scalar product with (zi\ in the last equation, we arrive at ^ = BY, where B = TT^ 
shows a circulant matrix structure (see Appendix [Bj. Using the diagonalization B = F^D^T^ 
of B, a Moore-Penrose pseudo-inverse can be computed as B + = ^^(D^ 1 )^ T^ N , and this allows 
to obtain f = B + ^> = T N {D- l fT ] N ^>. This last expression, by duality, can be interpreted as 
the convolution f = A * ^ between the data and the filter 



Remark 4.9. The relation between and T{k) is simply a "change of basis", but with non- 
orthogonal sets of generators {|^fe)} and {|5fc)}- Due to the particular choice of sampling points, 
the change of basis involves Fourier Transforms, and this can be interpreted as a convolution. 

Remark 4.10. For high spin values s » 1 (and N > 2s + 1), it is easy to realize that the filter 
(|36p acquires the simple form 



A(*)=!^(i+e-'/~ +0( i)) 



(37) 



where r = N — 2s. There is also a more manoeuvrable closed expression for the exact value of 
the filter zero mode 

^-i^sffl-'-wi:^! (38> 

n=0 n=0 
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where we have used the result of the Ref. [T7] concerning sums of the reciprocals of binomial 
coefficients. For large values of s one can also prove that 

2S£ (»)"' = »■ (39) 

71=0 

In the case of critical sampling all formulae are still valid, we only have to substitute N = 
2s + 1, the difference being that T is directly invertible and T _1 = 7^ + . The projector P is the 
identity, and A and B are both invertible. The reason for considering the case of oversampling 
is twofold: first, by its intrinsic interest leading to overcomplete frames, and second, in order 
to apply fast extensions (as FFT, see [I]) of the reconstruction algorithms it would be useful to 
consider N the smaller power of 2 greater or equal to 2s + 1. 

4.1.2 Under-sampling and critical sampling 

Let us suppose now that the number of sampling points is N < 2s + l. We shall see that, for N < 
2s + 1, we cannot reconstruct exactly an arbitrary function ip G 7i s but its orthogonal projection 
V> = Pn^ onto the subspace H s of Tt s spanned by the discrete set S = {{z^), k = 0, . . . , N} of 
CS. In other words, the restriction to this discrete subset implies a loss of information. 

This loss of information translates to the fact that the resolution operator A is not invertible 
and therefore we do not have a frame nor a resolution of the identity like in the previous 
subsection, see the discussion at the end of Sec. [2j But, since the set S is linearly independent, 
we can construct another operator, the overlapping kernel B = TT\ which is invertible and 
provides a partial reconstruction formula. In addition, the overlapping kernel operator has a 
circulant structure, and this provides a deep insight in the reconstruction process. 

Let us formalize again the previous assertions. 

Lemma 4.11. Let Q = {z^ = e 2mk ^ N , k = 0, . . . , N— 1} the discrete subset of the homogeneous 
space Q = SU(2)/U(1) = S 2 = C made of the N < 2s + 1 roots of unity. The pseudo-frame 
operator T : TL S — > given by T(ip) = {(zk\ip), Zk S Q} [remember the construction after 
Eq. (0|/7 is such that the overlapping kernel operator B = TT^ is an N x N hermitian positive 
definite invertible matrix, admitting the following eig en- decomposition: B = J-nDJ 7 ^, where 
D = diag(Ao, • • • , Ajv_i) is a diagonal matrix with = ^ Y2j=o ik+jN)' 9 being the ceiling of 
(2s + 1)/N. 

Proof. Let us see that the eigenvalues Afc of B = TT* are indeed all strictly positive and hence 
B is invertible. This can be done using RFM or taking advantage of the circulant structure of 
B (see Appendix [B]). With RFM we start with the expression of T = J~n,2s+i D 1 / 2 to obtain 
B = TT^ = Fn,2s+i D Ffr 2s+1 , which should be further worked on in order to fully diagonalize 
it. 

This can be done by using the "trick" mentioned in Appendix [A] consisting in enlarging the 
RFM J^7v,2s+i to ^/vm wnere M is the smaller multiple of N greater or equal to 2s + 1, and 
q = MJN is the ceiling of (2s + l)/iV (see Appendix [A}. In this way J-^m always contains an 
integer number of ordinary Fourier matrices J~n- 

Using this we obtain that B = J-j^ff J 7 ^^ where is the extension of D to a M x M 
matrix, and with a little of algebra the expression B = !Fn D J 7 ^ is obtained, where D = 
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diag(A , . . . , Ajv_i) and 

A * = E Xk + lN = tfTs E ik+m) ( 40 ) 

1=0 1=0 

All the eigenvalues are strictly positive and therefore B is invertible. ■ 
Following Sec. 2, we introduce the following: 

Lemma 4.12. Under the conditions of the previous Lemma, the set {\zk) = X^fcM) 1 (^ _1 )ifcl z z) > k = 
0, . . . , N — 1} constitutes a dual pseudo-frame for S, the operator P$ = T^T is an orthogonal 
projector onto the subspace of Jig spanned by S, where = T^B~ l is a (right) pseudo-inverse 
for T , and 

N-l N-l 

E \^}(z k \ = E \ z k)(zk\ = Ps (41) 

fc=0 k=0 

provides a resolution of the projector P$. 

Proof. If we define X r + = T* f B~ 1 it is easy to check that TT^ = In is the identity in C^. In 
the same way, Ps = T^~T is a projector since Pj = T+TT+T = T+T = P s and it is orthogonal 
pt = (Tt^- 1 T)t = rtB _1 T = Ps since B is self-adjoint. The resolution of the projector is 
provided by eq. (fT4"j). ■ 

Although the fully reconstruction of the original signal is not possible in the case of under- 
sampling, a partial reconstruction is still possible in the following sense. 

Theorem 4.13. (Partial reconstruction formula) Any function ip G TL S can be partially 
reconstructed from N < 2s + 1 of its samples (the data) ^(z^) = {zk\ip), at the sampling points 
Zk = e 2mk l N , k = 0, . . . ,N — 1, by the alias \iJj) = Ps\i>), by means of: 

N-l 

= (z$) = ^(zk)E(zz; 1 ), (42) 

k=0 

where 

N-l q-l 

= ^(i + **r - E v 1 E x p+inz p+in (43) 

p=0 1=0 

plays the role of a "sinc-type function". 

Proof. The proof follows the same lines as in Theorem l4.4i From the resolution of the projector 
(f4TI) . any ip G Tt s has an unique alias i/? = Ps^p which can be written as \ip) = ^( z k)\ z k)> 
and therefore ^(z) = (z\tp) = Y,k=o ®( z k)(z\zk)- Using that \z k ) = YjI=o^ \B~ l )ik\ z i) , we derive 
that 

N-l 2s s N-l q-l 

* 1=0 n=0 p=0 1=0 

(44) 

where eq. (j22|) and the orthogonality relations ([28]) have been used (but with N < 2s + 1, as in 
Appendix [B]). ■ 

Remark 4.14. As in the case of oversampling, eq. (|42p can be interpreted as a Lagrange- 
type interpolation formula, where the role of Lagrange polynomials are played by the functions 
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Lk(z) = ^(zz^ ), this time satisfying the proper orthogonality relations L^zi) = E(ziz^ ) = Sif,. 
The reason for this is that in the case of under-sampling there is not an overcomplete set of 
data, and therefore the "Lagrange functions" are orthogonal, although not complete. 
A partial reconstruction can also be obtained, in a natural way, from the "dual data" : 

Proposition 4.15. If we define the "dual data" as T(k) = (zk\ip), then they are related to the 
data *i>(k) = (zk\ip) through the convolution product: 

N-l 

r(jfe) = [A = A(fc -/)*(/), (45) 

1=0 

where A(/c) (the filter) turns out to be the discrete Fourier transform of 5 = (Aq 1 , . . . , A^,), 
where Xk are the eigenvalues of the overlapping kernel operator B; that is: 

N-l 

A(k) = [T N 6](k) = -= \^e~ i2 ™ k / N . (46) 

n=0 

Proof. The proof follows the same lines as in Proposition 14.81 with the difference that now 
B is not singular and there is no need for a pseudo-inverse. From = BT and using the 
diagonalization B = TjqDT^ of B, the inverse is directly B~ x = FnD' 1 ? 7 ^, and this allows to 
obtain f = B~ l $ = J^nFJ- 1 ^^. This last expression, by duality, can be interpreted as the 

convolution T = A * ^ between the data and the filter (l-KiD .B 
The comments made in Remark 14.91 also apply here. 

Again, a reconstruction in terms of the Fourier coefficients can be directly obtained by means 
of the (right) pseudo-inverse of the frame operator T: 

Corolary 4.16. The Fourier coefficients a m of the expansion \ip) = Y^m=-s®m\ s ' m ) °f the 
alias of any ip 6 H s in the angular momentum orthonormal basis B(H S ) can be determined in 
terms of the the data ^{z^) = (zk\ip) as: 

N-l N-l 

a n -s = § O V2 £ e 2 ^ N £ <B-%9{1) , n = 0, . . . , 2* . (47) 

k=0 1=0 

Proof. Taking the scalar product with (zk\ in the expression of we arrive at the system of 
equations 

2s 

Y,%na n -s = *(zfc), T kn = (z k \s,n - s), (48) 

n=0 

which can be solved by left multiplying it by the (right) pseudo-inverse of T, T r + = = T^B~ X . 
Using the expressions of B, given in Lemma 14.111 and the matrix elements 7fc n , given by the 
formula (|22p . we arrive at the desired result by noting that T^T = P$ and this acts as the 
identity on a n _ s . ■ 

Remark 4.17. Using vector notation this can be written as a = T^Ta = = T^B^ 1 ^, and 

this is even simpler in terms of the dual data, a = 7~tf = D 1 / 2 J 7 ^ 2s+l f . 

It is interesting to establish the connection between our results and others in the literature 
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Corolary 4.18. (Covariant interpolation) For < k < N — 1 define on Q the functions 
$fc(z) = (z\zk), z £ C. Let Co) • • • jCat-i be N complex numbers and Bki the overlapping kernel 
operator. Define on Q the function 



<&(*) 



. ZN-l] Co> • • • j Cn-1) z) 

( $i(z) 
Co £>0,0 



det(tf) 



det 



\ Gv-i £>tv- 



1.0 



Bo,N-l 

Bn-i,n-i ) 



(49) 



T/ien we have that: 

1. $(z) = (z\4>) for some <fi £ 7i s 

2. 3? is a solution of the interpolation problem, i.e.: ^(zf.) = Cfc> z = 0, . . . , N — 1. 

3. <3? is of minimal norm, in the sense that if $ is any other function on Q with 3?(z) = (z\4>), 
for some <ft £ 7i s , and &(z)-) = Ck> then \ \&\\ > H^H- 

Jf.. The interpolation procedure is invariant under left multiplication in G, in the sense that 
U(g)BU( g y = B and 

${gz , gzN-r, Co, • • • , Cn-wqz) = $(zq, . . . , z N -r, Co, • • • , Ov-i; z), 

(gz denotes the natural action of the group G on its homogeneous space Q = G/H) so 
that the left- displaced interpolation problem ^(gz^) = Cfc is solved by the function $(z) = 



Proof. It is a direct consequence of Theorem 14.131 if we identify the data Cfe = ^>(zk), &(z) = 
&(z) and Ylk=o (^^kl^liz) = ^{zz^ 1 ). The fact that $ is of minimal norm is a direct conse- 
quence of the orthogonality of the projector P$. The invariance under left multiplication is a 
consequence of the invariance of the overlapping kernel B under left multiplication. ■ 



4.2 Several spin case 

The case of several spins, i.e. band limited functions, is more involved than the single spin case, 
and it is not so easy to select the sampling points in such a way that an explicit expression for 
the inverse of the resolution or overlapping kernel operators be available. 

Denoting bj| H {J) = ® J S=0 H S the Hilbert space of band-limited functions, up to spin J, the 
set of coherent states can be defined in an analogous way to the single spin case. 

First, let us denote by U^(z,z) = (& J S=0 U S (z, z) the unitary and reducible representation 
of SU(2) acting on TL^ J \ where U s (z,z) stands for the unitary and irreducible representation of 
spin s. The Hilbert space TiS J ' has an orthogonal basis given by {|s,m)}, in such a way that 
^H( J) = jxx S^=o Yltn=-s I s ' m )( s i m \ 1S a resolution of the identity. Selecting the fiducial vector 
1 7) J = ©g=o I s ) s >, the set of coherent states is defined as \z) J = U^ J \z, z)\i) J ■ 

§ For the time being we shall restrict ourselves to integer values of spin, in order to compare with standard 
Fourier Analysis on the sphere. 
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The CS overlap, for the several spins case, is now: 

s=0 ' 

The first, naive choice, of sampling points would be the N roots of unity, z^ = e 2lT%k / N , where 
now N = dinxH^ = ( J + l) 2 in order to have critical sampling. In this way the operators A^ J ' 
and B^ J ' would have nice structure and their inverse matrices would be easily computed. 

However, the following negative result prevents us from proceeding in this way: 

Proposition 4.19. For N > 2J + 1, the overlapping kernel operator B^ J ' has range 2J + 1. 
Proof. Let A s , T s , B s and D2 S +i the eigenvalues, frame, overlapping kernel operators and 
diagonal matrix appearing in the previous sections corresponding to angular momentum s. 
Then the frame operator T( J > ; 7i^ J ' — > can be written as a N x ( J + l) 2 matrix given by 

■^oL) = T k s n = Fn,2s+i vVs+V 

Then £( J ) = T^T^ = B +B 1 +- • -+B J = T N (d\ +d\+- ■ -+d\ j+1 ) = F N d\ j+1 
where -D2J+1 is a diagonal matrix with eigenvalues: 

1 J 

~ Xn = JTT S A "' n = 0,l,...2J (51) 

s=(n-l)/2 

where (n — l)/2 stands for the ceiling of ^^i^-B 

The proof could also have been done using the circulant structure of B^ J \ which can be 
written as Bff = Ci- n , where now 

s=0 

and computing its eigenvalues as in Appendix iBl 

Therefore, putting all sampling points in the equator of the Riemann sphere is not a good 
choice, and other alternatives should be looked for. The problem is that other choices of sampling 
points lead to resolution operators with less structure and therefore without the possibility of 
having an explicit inverse. 

Another possibility is to use an equiangular grid in (9,4>), as the one used in [3J. If 
(Oj, <pk) = jfk) , j,k = 0,1, ... ,N — 1 is a grid of A^ 2 points in the sphere, where A^ > 
J + 1, the corresponding points in the complex plane by stereographic projection are given by 
z k = e %< ^ k tan = e l ~ k t&n(^j) = rje l ~ k . However, it can be checked that in this case the 
resolution operator is also singular. 

We shall follow a mixture of both approaches, consisting in using as sampling points the 
2s + 1 roots of (r s ) 2s+1 , for s = 0, 1, . . . , J. Here r s is a positive number depending on s, in such 
a way that if s 7^ s' then r s ^ r s >. Thus, we shall continue to use N = (J + l) 2 sampling points 
but distributed in circles of different radius. In the Riemann sphere, these would be distributed 
in different parallels, one for each value of spin. These points are given by: 

z^j = r s e 2a + 1 , s = 0,..., J, m = 0,...,2s, 
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where s denotes spin index and m the third component of spin. 

The overlapping kernel operator for this choice of sampling points is given by: 



1 J 

B (a,b) _ /Ja)\ (b)\J_ 1 V 
^m,n \ m \ n I j . j / j 

s=0 



(l+rl)(l+rl) 



r, ■ n(2a + l)-m(2b+l) \ 2 \ S 

1 + r a r b e [25+rjps+rj 



1 if z (a) - 



1 m a "h , otherwise 



J+i 



1 WjTJ..' 



(52) 



where 



ab _ (1 + r a r b e 2 ™ Wd^+i, ' )2 



.re(2a + l)-m(2i)+l) 



(l+^)(l+r b 2 ) 



is the multiplier of a geometric sum. 



The overlapping kernel operator Bm',n is an hermitian matrix having the following structure: 



(53) 



circ(l) 




S02 




Bofc 


-Bo fc+i 


... \ 


Bl 


circ(CW,CW,CW) 






Bifc 


Bi fc+i 




4i 




circ(C ( l ) 2 ^...,C| 2, ) 




i?2fc 


S2 fc+i 


















4b 


4* 


4, 






Sfc fc+i 




^0 fc+l 




Rt 

D 2 fc+l 




B\ fc+i 





















M (1 + r 2 e 27rin/(2s+l))2 

The diagonal blocks are circulant matrices of dimension 2s+l, with Cn = 1 ^tt> 

(1 +r^) 2 

and the non-diagonal blocks -E> p(? are matrices of dimension (2p + 1) x (2q + 1). 

The overlapping kernel operator S^'n is not a circulant matrix, not even a block circulant 
matrix, therefore the computation of its inverse needed for the reconstruction formula must be 
done numerically. Even the checking that it is non-singular must be done numerically. 

For different choices of r s , we have checked that the resolution operator is invertible, and 
therefore the reconstruction formula can be used, although nice expressions like (|30p or (|42p are 
not available. 

An interesting feature observed for all choices of r s is that the overlapping kernel operator 
does not seem to approach the identity matrix when J — * 00, as one would naively expect, i.e. 
in the limit of large J, the frame does not approach to a tight frame. Even more, it seems that 
the overlapping kernel operator B^'n approaches to a singular operator since its eigenvalues 
tend to zero as J grows. This is not a serious drawback, however, since we are interested only 
in band-limited functions on the sphere and therefore J remains always finite. 
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5 Connection with the Euler angle picture 



We have provided reconstruction formulas for Majorana functions ^f(z) from N of its samples 
= {Z)M>) a t the sampling points the Riemann sphere. The advantage of using 

this "complex holomorphic picture" , instead of the standard "Euler angle picture" , is twofold: 
firstly we can take advantage of the either diagonal or circulant structure of resolution and 
overlapping kernel operators, respectively, to provide explicit inversion formulas and, secondly, 
we can extend the sampling procedure to half-integer angular momenta s, which could be useful 
when studying, for example, discrete frames for coherent states of spinning particles in Quantum 
Mechanics. 

Moreover, for integer angular momenta s = j, we could always pass from one picture to 
another through the Bargmann transform (|25|) . Indeed, let us work for simplicity in the critical 
case N = 2j + 1 and let us denote by <£fc = {9q, (pkl^) the samples of the function ([20]) . in the 
Euler angle characterization, at the sampling points 9q ^ 0, it and = — k = 0, . . . , N — 1 
(i.e., a uniformly distributed set of ./V points in a parallel of the sphere S 2 , but counted clockwise). 
Denoting by: 

K kl = K(9 ,cj )k ;z l ) = (e Ak\zi) = ^M e ^^ s i nJ( o )(i + 2cot(0 o )e^( i - fc )-e^('- fc ))^ (54) 

a discrete N x N matrix version of the Bargmann transform (|25p . and inserting the resolution 
of the identity (|26p in (#o>0jfc|V'}> we easily arrive to the following expression: 

AT-l 

$k = ^iB^m, (55) 

l,m=0 

which relates data between both characterizations or pictures through the CS transform and CS 
overlap matrices K, and B in ([53]) and (fT3j) . respectively. 

Except for some values of 9q ( see later in this section) , the transformation (|55p is invertible 
and explicit formulas of /C _1 are available. Actually, K, can be written as the product K, = AQ, 

A *p = ^Pf e'^W'CWfcp, Q P i = (1 +2cot(0 o )e^«-*>) - e *$(*~P))> = qi _ p , ( 56 ) 

of a diagonal matrix A times a circulant matrix Q, which can be easily inverted (following the 
procedure of Appendix[B|) as = Fn^ 1 ^]^, where = diag(a;o, • • • , wjv-i), with eigenvalues 

JV-l j v 

"fc = £ 1ne- l % kn = 7V^^ / (-l) r ©(2cot(0o)r r , (57) 

n=0 p=0 r=0 

where the prime over ]TJ implies the restriction p + r = k. Therefore, we can also obtain data in 
the holomorphic characterization, from data in the Euler angle characterization, through 
the formula: 

* = BQ^A- 1 ® = TnDQ-^^A-^ (58) 

which can be seen as a convolution ^ = O * $' of the re-scaled data $' = A~ 1( £ times the filter 
= FnQ, with 9 k = Afc/cjfc the quotient of eigenvalues of B and Q. 
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Note that there are values of 8q f° r which fC is not invertible. Such is the case of Oq = 7r /2 
(the equator), for which u)k = N(—l) k / 2 ( k J / 2 ) if k even, and zero otherwise. Let us show that 
this situation is linked to the fact that general functions (|20p in the Euler angle picture can not 
be reconstructed from its samples on a uniformly distributed set of N points in the equator 
of the sphere. Indeed, let us insert this time the resolution of unity hj+i = Yjm=-j \j) m }(j> m \ 
in (0 o ,0a#), witn W) = YHn=-j a m\j, m ), which results in: 



E J^jYf\0 o ,cP k )a m , (59) 

m=-j 

where we have used the definition ffTU]) . Denoting J4n(#o) = \f$Yj l ~'' (@o, 4>k) and knowing 
from Remark 14.71 that the Fourier coefficients a n —j are given in terms of data VP^ trough a = 
we arrive to a variant of the formula (|55p : 

$ = y(9 )D-y 2 ^, (60) 

which again connects data between both pictures. Knowing that spherical harmonics can be 
expressed in terms of associated Legendre functions Pp by: 

Yp{e,<f>) = e iw ^Pf(cos0), (61) 
whose value at the equator 6q = ir/2 is given in terms of Gamma functions as: 

P m = 7=<*» (Uj + rn)) T r { ,\ J + \ m+ J ) y (62) 

we immediately realize that y(n/2)kn = for n odd. In other words, for #o = tt/ 2 , the 
reconstruction process in the Euler angle picture fails unless we restrict to the subspace of 
functions ip with null odd Fourier coefficients (i.e., a n _j = for n odd). 



A Rectangular Fourier Matrices 

Let N, M G N, and let T NM be the N x M matrix: 

(F N M)nm = -^=e i2nnm / N , n = 0, 1, . . . , JV - 1 , m = 0, 1, . . . , M - 1 . (63) 
v N 

We shall denote these matrices Rectangular Fourier Matrices (RFM). For N = M we recover 
the standard Fourier matrix Tn- Let us study the properties of these matrices in the other two 
cases, N > M and N < M. 



A.l Case N > M 

The case N > M is the one corresponding to oversamplig, and it is the easiest one, since it is 
very similar to the N = M case. Let us first introduce some notation. 
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Let lnm '■ C 



M 



be the inclusion into the first M rows of C (i.e. padding a M 



vector with zeros). And let pmn '■ C 
truncating a iV- vector). 

It can be easily checked that: 



N 



C the projection onto the first M rows of C (i.e. 



Pmn lnm = Im , i>nm Pmn = Pm : - 

where Im stands for the identity matrix in C M . 

Given the square matrices A and B acting on and C M , 
matrices A^ and through the commutative diagrams 



Im 












(64) 



NxN 



respectively, we define the square 



C 



ill 



A i 



(■JVM 



c 



N 



c 



M 



c 



Pmn 



N 



c n 

Pmn 
(TM 



* C N 

l>NM 

C M 



(65) 



The matrix = Pmn Alnm is the truncation of A to a M x M matrix and the matrix 



B^ = lnm B pmn = 



B 












is the padded version of B. Also note that Pm = (I. 



m 



it 



NxN 



From these definitions the following properties of the Rectangular Fourier Matrices are de- 
rived: 



•F NM = Fn l^NM , F 



•t _ 



NM 



PMN F 



t 



N 



Fmm Fnm 



Im , 



Fnm F n m = Fn Pm F n . 



(66) 



A. 2 Case N < M 



The case N < M is the one corresponding to undersamplig, and it is not as easy as the N > M 
case. The same definitions as in the previous case also apply here, but interchanging the roles 
of and M. Thus, in this case we have: 



Pnm imn = In , t-MN Pnm = Pn = 
The RFM in this case now read: 



In 












= (In)' 



(67) 



MxM 



Fnm 



{ F N I F N | i H mcs | Fn \ F Np ) 
/ F N \ 



F 



NM 



N 



q times 



F, 



r 



N 



(68) 



where p = M mod N and q = M div N. 
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Instead of working with these matrices, it is more convenient to "complete" them so as to 
have an integer multiple of Fourier matrices. Le M be the smaller multiple of N greater or equal 
to M, and q = M/N the ceiling of (2s + 1)/N. Note that 



q p = 
9 + 1 P^O 



Then 



NM 



T N \T N \ 9 *P»* | T N ) 



(69) 
(70) 



a similar expression for J~ N ^ is obtained, and 



NM l MM i F NM — PmM -^NM 

qIn + Fn Pp F N ' 3~nm Fnm = (Im)^ = hi 



rt 



Fnm F\ 



t 

WM 



where 



In 


In 


q times 


In \ 




In 


In 




In 




'. q times 










In 


In 




In ) 


MxM 



(71) 



(72) 



B Circulant Matrices 

The overlapping kernel operator B has a circulant matrix structure which gives a deep insight 
into the process taking place and we may take advantage of this fact to diagonalize it in the 
case of undersampling where RFM are more difficult to handle. 
Indeed, note that: 



B k i = (z k \zi) = ^ (l + e 2 ^- fc )/^) 23 = d_ k , k,l = 0,...,N-l. 
where C n = ^ (l + e 2mn / N ') 2s ; which shows a circulant matrix structure 



(73) 



where 



B = circ(C ,Ci, . . .,C N -i) 



/ 1 ... \ 



/ C Ci 

C-N-l Co 



C 2 



Cn-i \ 

Cn-2 

Co J 



N-l 

^ c.if = p c (n), 

J=0 



(74) 



n 



(n" = /at, 

i 

V 1 ••• °/ 

is the generating matrix of the circulant matrices and P c (t) is the representative polynomial of 
the circulant (we put II = In). According to the general theory (see e.g.|12j). every circulant 
matrix is diagonalizable, whose eigenvectors are the columns of the Vandermonde matrix Vjv = 



n* = n f = ir 



n^ 1 ), 
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polynomial a: 



V(zq, . . . , zn-j) = vNJ-^ and whose eigenvalues can be computed through its representative 

Zk) = 
5-1 



N-l N-l 2s 

\ k = Pc (z k ) = ClZk 1 = 2 ~ (2S) E E ( 2 ;)e M(2s -" )/J, e- 2nM / w 

1=0 1=0 n=0 



= £s-Y,ik+w),k = Q,...,N-l, (75) 

1=0 

where q is the ceiling of (2s + 1)/N and we have used the orthogonality relation (|28p . although 
in this case, since N < 2s + 1 there can be more terms in the sum. All of them are strictly 
positive, and it is easy to proof that B = F^DF^, where D = diag(Ao, • • • , Aat_i). ■ 
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